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ABSTRACT 

One of the major challenges to identification of the 3.3, 6.2, 7.7, 8.6, and 11.3/um 
interstellar infrared (IR) emission bands with polycyclic aromatic hydrocarbon (PAH) 
molecules has been the recent detection of these bands in regions with little ultraviolet 
(UV) illumination, since small, neutral PAH molecules have little or no absorption at 
visible wavelengths and therefore require UV photons for excitation. 

We show here that our "astronomical" PAH model, incorporating the experimental 
result that the visual absorption edge shifts to longer wavelength upon ionization and/or 
as the PAH size increases, can closely reproduce the observed infrared emission bands 
of vdB 133, a UV-poor reflection nebula. 

It is also shown that single-photon heating of "astronomical" PAHs in reflection 
nebulae near stars as cool as T e g- = 3000 K can result in observable emission at 6.2, 7.7, 
8.6, and 11.3/um. Illustrative mid-IR emission spectra are also calculated for reflection 
nebulae illuminated by cool stars with T e g = 3500, 4500, 5000 K. These will allow 
comparison with future Space Infrared Telescope Facility (SIRTF) observations of vdB 
135 (T eS = 3600 K), vdB 47 (T cff = 4500 K), and vdB 101 (T eff = 5000 K). 

It is also shown that the dependence of the 12/im IRAS emission relative to the total 
far-IR emission on the effective temperature of the exciting star is consistent with the 
PAH model expectation for 3000 K < T eff < 30000 K. 

Subject headings: dust, extinction — infrared: ISM: lines and bands — ISM: individual 
(vdB 47, vdB 101, vdB 133, vdB 135) — reflection nebulae 



1. Introduction 

Since their first detection in the planetary nebulae NGC 7027 and BD+30°3639 (Gillett, For 
rest, &: Merrill 1973), the distinctive set of infrared (IR) emission features at 3.3, 6.2, 7.7, 8.6 
and 11.3/xm have been observed in a variety of objects with a wide range of physical conditions 



2 



including planetary nebulae, protoplanetary nebulae, reflection nebulae, HII and ultracompact HII 
regions, circumstellar envelopes, the diffuse interstellar medium (ISM) of the Milky Way Galaxy, 
and external galaxies (see Tielens et al. 1999 and Sellgren 2001 for recent reviews). Although the ex- 
act nature of their carriers remains unidentified - they remain known as "the Unidentified Infrared 
(UIR) bands" - it is now widely thought that they originate from some sorts of aromatic hydro- 
carbons. Various carbonaceous materials have been proposed as the UIR band carriers. In general, 
they can be categorized into two classes: (1) pure, free-flying aromatic molecules - polycyclic aro- 
matic hydrocarbon molecules (PAHs; Leger & Puget 1984; Allamandola, Tielens, & Barker 1985); 
(2) carbonaceous grains with a partly aromatic character - hydrogenated amorphous carbon (HAC; 
Duley & Williams 1981; Jones, Duley, & Williams 1990), quenched carbonaceous composite (QCC; 
Sakata et al. 1990), coal (Papoular et al. 1993), fullerenes (Webster 1993), and surface-graphitized 
nanodiamonds (Jones k, d'Hendecourt 2000). 

In most current models, the UIR band emission involves three sequential steps: (1) excitation 
by absorption of an energetic starlight photon (usually ultraviolet [UV]); (2) rapid (~ 10~ 12 — 
10~ 10 s) redistribution of all or part of the absorbed photon energy over all available vibrational 
modes; (3) radiative relaxation via IR fluorescence. 

Among the existing proposed carriers, the PAH model is gaining increasing acceptance because 
of (1) the close resemblance of the UIR spectra (frequencies and relative intensities) to the vibra- 
tional spectra of PAH molecules (e.g. see Allamandola, Hudgins, & Sandford 1999) and (2) the 
ability of a PAH molecule to emit efficiently in the UIR wavelength range following single photon 
heating (Greenberg 1968; Leger & Puget 1984; Allamandola, Tielens, & Barker 1985; Draine & 
Li 2001). In contrast, larger carbonaceous grains are unlikely to be efficient UIR emitters since 
the timescale for the absorbed photon energy to diffuse (~ 10~ 9 s) is much shorter than the IR 
vibrational emission timescale (~ 0.1s; see Tielens et al. 1999). 

Early observations of the UIR emission bands were made in regions with strong UV irradiation, 
and PAH excitation was expected since all PAH species are strongly absorbing in the vacuum 
ultraviolet (A < 3000A). The PAH UV excitation model has recently been challenged by the ISO 
{Infrared Space Observatory) detection of the UIR bands in vdB 133, a reflection nebula illuminated 
by a binary system with little UV radiation (see Uchida, Sellgren, & Werner 1998). The UIR 
spectrum of this UV-poor region closely resembles those observed in sources with much harsher 
UV environments (Uchida et al. 2000). This appears to be in conflict with the view that PAHs are 
primarily excited by UV photons, as would be expected based on laboratory studies showing that 
the absorption by small, neutral PAHs has a sharp cut-off in the UV, with little or no absorption 
in the visible (see Sellgren 2001 for a review). Uchida et al. (2000) obtained 5-15/mi ISOCAM 
spectra of 6 reflection nebulae and found no systematic spectroscopic differences despite values of 
T eff ranging from 6800 K to 19000 K. 

The PAH electronic absorption edge is known to shift to longer wavelength with increasing size 
and/or upon ionization (see Allamandola et al. 1989, Salama et al. 1996 and references therein). 
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While the largest experimentally studied PAH molecule to date is dicoronylene C48H20 (Allaman- 
dola, Hudgins, & Sandford 1999), astronomical PAHs are believed to be larger (e.g., the mean size 
for the Milky Way PAHs is ~ 6A, corresponding to Nq ps 100 [Li & Draine 2001a]). The astro- 
nomical PAH model - with the size/ionization dependence of the PAH absorption edge taken into 
account (see §A2 in Li & Draine 2001a) - is successful in explaining the observed mid-IR spectra 
of the Milky Way diffuse ISM (Li & Draine 2001a), the quiescent molecular cloud SMC Bl#l 
in the Small Magellanic Cloud (Li & Draine 2001b), and the UIR band ratios for a wide range 
of environments ranging from reflection nebulae, HII regions, photodissociation regions (PDRs), 
molecular clouds in the Milky Galaxy to normal galaxies, starburst galaxies, and a Seyfert 2 galaxy 
(Draine & Li 2001). 

In this Letter we show that the astronomical PAH model is consistent with the observed 
UIR emission from UV-poor environments. In §2 we verify that the astronomical PAH model 
can quantitatively reproduce the vdB 133 UIR spectrum. We further demonstrate in §3 that 
the UIR bands are also expected for reflection nebulae which are even more UV-poor than vdB 
133, and we provide model spectra for comparison with future Space Infrared Telescope Facility 
(SIRTF) observations of reflection nebulae near cool stars. In §4 we show that the predicted T e g 
dependence of the ratio of the IRAS 12/j.m emission to the total far-IR emission is also consistent 
with observations. Our results are discussed in §5, and our conclusions are summarized in §6. 



2. vdB 133 

The reflection nebula vdB 133, located in the Cyg OBI association at a distance al ps 1500 pc, 
is illuminated by a binary system consisting of HD 195593A (spectral type F5Iab; effective temper- 
ature T cS ps 6800 K; absolute magnitude My = -6.6 ± 0.4) and HD 195593B (spectral type B7II; 
T cff Ps 12000 K; My = -4.0 ±0.4) (see Uchida et al. 2000 and references therein). Following Uchida 
et al. (1998), we approximate the vdB 133 radiation field by a linear combination of two Kurucz 
model atmosphere spectra (Kurucz 1979). The starlight intensity incident on the reflection nebula 
is assumed to have Uq ps 135, where Uq is the ratio of the 912A-l^m energy density relative to the 
value for the MMP interstellar radiation field (4.80 x 10~ 13 erg cm -3 ). 

Let f\ be the fraction of the nonionizing stellar luminosity L\ radiated at wavelength < A: 

f A ° LydX' 

f — J912A A (1\ 

JX - r o LydX' 1 ' 

J 912A A 

This binary system, 1 with /o.2^ m ~ 1-7%, has the lowest value of /o.2^m for reflection nebulae with 



1 HD 195593B, with / . 2Mm « 32%, contributes « 80% of the total flux shortward of 0.2^m; HD 195593A, with 
/o.2 M m ~ 0.4%, contributes w 95% of the total flux shortward of 1/ttm. 



4 



detected UIR emission bands. 2 Despite this, the vdB 133 UIR features are very similar to those 
observed in regions with much harder radiation fields (Uchida et al. 2000). 

We model the vdB 133 UIR features as emitted from a mixture of neutral and ionized PAH 
molecules. We use the PAH absorption cross sections estimated by Li & Draine (2001a) which take 
into account the size/charge dependence of the PAH long wavelength absorption edge (see Figure 
1). We adopt a log-normal size distribution for the PAHs, characterized by three parameters: do, 
a, and 6c; and o respectively determine the peak location and the width of the log- normal 
distribution, and 6c is the total amount of C atoms (relative to total H) locked up in PAHs. As 
always, the term "PAH size" refers to the radius a of a hypothetical spherical grain with the same 
carbon density as graphite (2.24 gem -3 ) and containing the same number of carbon atoms Ac: 
a = 1.288A C /3 A. 
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Fig. 1. — Absorption cross sections (per C atom) estimated for astronomical PAH molecules containing iVc = 
20, 35, 50, 70, 100, 200 C atoms, for both charged (upper panel) and neutral (lower panel) states. 

The PAH ionization fraction depends on the radiation intensity, the electron density n e , and the 
gas temperature T gas . Assuming T gas = 100 K, we estimate the PAH ionization fractions in vdB 133 
using the rates for photoelectric emission and electron capture recently discussed by Weingartner 



2 For comparison, NGC 2023 is illuminated by HD 37903, with T e g = 22000 K and / .2 M m ~ 68%; the average 
interstellar starlight background (Mathis, Mezger, & Panagia 1983) has /o.2 M m = 4.7%. 
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k Draine (2001). 

We employ the "thermal-discrete" method (Draine &; Li 2001) to calculate the temperature 
distribution functions for PAHs subject to single-photon heating. Our model assumes a single 
radiation field: attenuation of the illuminating starlight in the dust layer is neglected. If the dust 
layer is optically-thick to the exciting radiation, we may approximate this by an optically-thin layer 
with a column density Ah such that the optical depth to the exciting radiation would be unity. For 
a standard dust mixture and A = 4000 A (say), this would correspond to Ah ~ 1-3 x 10 21 cm -2 . 

Figure 2 plots the theoretical spectra calculated from models with ao = 3.5A, a = 0.4, and 
abundances (A H /10 21 cm" 2 )& c = 38,33,29ppm for U /n e = 250,500,1000 cm 3 respectively. We 
see that the Uo/n e = 500 cm 3 model (n e = 0.3 cm -3 ) closely reproduces the observed UIR spectrum 
using the standard Milky Way PAH size distribution (ao = 3.5A, a = 0.4), a realistic radiation field, 
and PAH abundances (6c ~ 33ppm if Ah = 10 21 cm -2 ) comparable to the Milky Way abundances 
(b c ~ 45ppm). 




A (/am) 

Fig. 2. — Predicted IR emission spectra for PAHs in vdB 133 in comparison with the ISOCAM observational 
spectrum (filled circles; Uchida et al. 1998). The PAH mixture is modeled as a log-normal distribution of ao = 3.5A, 
a = 0.4, b c » 38ppm x (10 21 cm _2 /Ar H ) for U /n e = 250cm 3 (dotted line), b c « 33ppm x (10 21 cm" 2 /JV H ) for 
U /n e = 500cm 3 (solid line), and b c » 29ppm x (10 21 cm^/Nn) for U /n e = 1000cm 3 (dashed line). We take the 
observed UIR band widths (see Table 7 of Li & Draine 2001). 
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3. Cool Reflection Nebulae 

We have seen in §2 that the astronomical PAH model is successful in explaining the UIR 
spectrum observed in vdB 133. In this section we will show that reflection nebulae illuminated by 
cooler stars are also capable of exciting PAHs sufficiently to emit at the 6.2, 7.7, 8.6 and 11.3/xm 
UIR bands. For illustrative purpose, we consider reflection nebulae illuminated by stars with 
T e ff = 3000, 3500, 4500, 5000 K. The latter 3 cases will serve as a direct comparison basis for SIRTF 
observations of vdB 135 (T eff = 3600 K), vdB 47 (T cff = 4500 K), and vdB 101 (T cff = 5000 K) 
(Houck 2001). 3 

Lacking prior knowledge of the PAH size distribution and the physical environmental parame- 
ters such as n e , T gas , and the starlight intensity Uq, we adopt the Milky Way PAH size distribution 
(i.e., a log-normal distribution with oq = 3.5 A and a = 0.4 [Li & Draine 2001]) and we take 
T g3S = 100 K. We take Uq = 100 as a representative number since 30 < Uq < 300 for typical 
reflection nebulae (see Sellgren et al. 1990); the calculated UIR emission is essentially proportional 
to Uq for Uq < 10 4 (see Figure 3a). In all models we take the number of carbon atoms relative to 
hydrogen in PAHs to be be = 45ppm, the same as in the Milky Way diffuse ISM (Li & Draine 
2001). 

For stars with T c g > 3500 K we take the Kurucz atmospheric model spectra to approximate 
the radiation fields. Since the Kurucz model is not available for stars cooler than 3500 K, we take 
the model spectrum of Allard et al. (1997) for T cS = 3000 K. 

In Figure 3a we plot the mid-IR spectra calculated for PAHs heated by a T c g = 3000 K star 
with Uq = 1,100, 10 3 ,10 4 and U /n e = 10 4 . The predicted PAH IR bands are prominent, with 
a much stronger 11.3/Um C-H out-of-plane band in comparison with those in regions illuminated 
by much hotter stars (e.g., NGC 2023 [T cff = 22000 K], Orion Bar [T cS = 37000 K], M17 [T cff = 
45000 K]; Verstraete et al. 2001). This is due to the lack of UV photons capable of causing 
photoelectric ejection of electrons and thus PAHs in such a cool region mainly remain neutral even 
for UQ/n e = 10 4 . As expected from the single-photon heating scenario, the PAH IR features are 
almost indistinguishable for Uq < 10 4 except the Uq = 10 4 model has a rising wing for A > 10/iim. 
This is because for Uq = 10 4 the > 20 A dust grains already attain an equilibrium temperature 
T » 100 K. 

In Figure 3b,c,d we plot the model spectra for Uq = 100, n e = 0.01, 100 cm -3 , and T e g = 
3500, 4500, 5000 K, respectively. It is seen that all spectral shapes are almost the same except the 
absolute emissivity level (at fixed Uq) increases with increasing of T e g. This is because hotter stars 
radiate a larger fraction of their power at short wavelengths where PAHs are strongly absorbing. It 
is also seen from the enhanced 11.3/xm band relative to the 7.7fim C-C stretching band that a large 



3 Attempts by Uchida et al. (2000) to search for UIR bands in vdB 101 and vdB 135 using ISOCAM spectroscopy 
were unsuccessful due to the faintness of these nebulae. 
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Fig. 3. — Predicted 5-14/im emission spectra for reflection nebulae illuminated by cool stars with an effective 
temperature of T cff = 3000 K and starlight intensity U = 1, 100, 1000, 10 4 (a), T e g = 3500 K and U = 100 (b), 
T cff = 4500 K and U Q = 100 (c), and T e g = 5000 K and U = 100 (d). In all models the PAH size distributions 
are taken to be that of the Milky Way diffuse ISM: a log-normal size distribution with ao = 3.5A, a — 0.4, and 
be = 45ppm. We take n e /n H = 1.2 x 10" 4 , T gas = 100 K, and U /n e = 1.0 cm 3 (solid lines) or Uo/n e = 10 4 cm 3 
(dotted lines) for calculating the PAH charging. 



fraction of PAHs in those regions are neutral. 4 This is the major difference between the UIR features 
excited by cool stars with those excited by hot stars. 5 It is hoped that SIRTF observations of these 
reflection nebulae will provide both spectra of the UIR bands and a full spectral energy distributions 
to allow the illuminating intensity Uq to be well-determined, thereby allowing comparison with the 
model spectra in Figure 3. 

In Table 1 we give the PAH mid-IR emission convolved with the 5.8/um and 8.0/iin SIRTF 
Infrared Array Camera (IRAC) bands calculated for the Milky Way PAH mixture excited by cool 



4 The probability for a small PAH molecule being neutral is larger for Uo/n e = 10 4 cm 3 than that for Uo/n e = 
1.0 cm 3 . This is because in regions devoid of UV photons the charging of a neutral PAH is dominated by electron 
capture and thus PAHs are more likely to be negatively charged when Uo/n e is small. 

5 The vdB 133 nebula is not lack of photoelectron-ejecting photons. The hot companion HD 195593B (T c g = 
12000 K) of the exciting binary system contributes » 80% of the total 6-13.6 eV photons. This is why the vdB 133 
UIR spectrum resembles those excited by hot stars. 
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stars with T cS = 3000, 3500, 4500, 5000 K. Figure 3 and Table 1 may provide guidance for the 
coming SIRTF observation of cool reflection nebulae (Houck 2001). We stress, however, that the 
actual mid-IR spectra may differ from those in Figure 3 and Table 1 since the PAH size distribution, 
the PAH abundance, and the environmental parameters n e , T gas and Uq may differ from those 
adopted here. 

Table 1: PAH emissivity averaged over the 5.8/xm and 8.0/jm SIRTF IRAC bands for the Milky 
Way mixture illuminated by cool stars for Uo/n e = 1.0 (10 4 )cm 3 . 
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2.66 (2.01) 
10.1 (7.05) 


4.05 (3.17) 
14.7 (10.5) 



4. The 12/im Emission vs. Radiation Hardness 

To study the carrier of the 12[im emission and its excitation mechanism, Sellgren et al. (1990) 
investigated how the ratio [A/A(12/um)]iRAs/^boi(FIR) depended on the effective temperature T e g 
of the exciting stars for 24 reflection nebulae, where [A/A(12//m)]rRAS is estimated from the IRAS 
12/xm filter, and the total far-IR surface brightness Iboi(FIR) was estimated by fitting a blackbody 
multiplied by a 1/A emissivity to the IRAS 60 and lOO^rn surface brightnesses. They found that 
A7 A (12/xm)/J bo i(FIR) is basically independent of T cS for 5000 K < T cS < 33000 K. If PAHs are 
the 12/iin carrier and excited only by UV photons, one would expect [A/A(12/im)]iRAs/-^boi(FIR) 
to show a sharp drop for cool stars (T e g < 10000 K) which are relatively devoid of UV radiation 
(Sellgren et al. 1990). Therefore, one can conclude that either PAHs are not responsible for the 
12/xm emission or PAHs are excited by photons of a wide range of wavelengths. Since it has been 
well established that the UIR features are the dominant contributor to the 12//m emission, the 
former hypothesis can be ruled out. Therefore, it must be possible for photons longward of the UV 
to excite the UIR carrier. Is this consistent with our understanding of PAHs? 

To test the astronomical PAH model, we calculate [A/>(12/im)]TRAs/-^boi(FIR) f° r reflection 
nebulae illuminated by stars with T cfr =3000, 3500, 4500, 6000, 8000, 10000, 15000, 22000 and 
30000K and starlight intensity Uq = 1, 100, 1000, assuming Uo/n e = 100cm 3 and using the Milky 
Way dust mixture. As shown in Figure 4a the model results are in good agreement with obser- 
vational data except four nebulae: vdB 135 (T cS = 3600 K), vdB 42 (T cS = 4200 K), vdB 47 
(T e fj = 4500 K), vdB 35 (T c g = 4900 K). It is worth noting that these nebulae are the four faintest 
among Sellgren et al.'s sample of 24. It is possible that Iboi(FIR) was underestimated for these cool, 
faint nebulae as a consequence of the method of obtaining Jboi(FIR) because the 60 and lOO^m 
emission preferentially account the warm grain component while the bulk of the far-IR emission is 
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probably longward of 100/im for cool reflection nebulae. 

A better way to characterize the fraction of total IR emission emitted at 12/nm is to replace 
Iboi(FIR) by iboi(tot) = J °° I^dA. In Figure 4b we show the model predicted [A/^(12/mi)]rRAs/^boi(tot) 
as a function of T e g . 
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Fig. 4. — Ratio of the 12^m emission to the total far-IR emission estimated from the 60 and 100/^m IRAS fluxes 
(upper panel; see text) and to the actual total absorbed energy (lower panel). Observational data (filled circles) 
are taken from Sellgren, Luan, & Werner (1990). Model results are for reflection nebulae illuminated by stars with 
T cff = 3000, 3500, 4500, 6000, 8000, 10000, 15000, 22000, 30000 K and starlight intensity Uo = 1, 100, 1000 times that of 
the 912A— ljLtm MMP interstellar radiation field, assuming Uo/n e = 100 cm 3 and using the Milky Way dust mixture. 



5. Discussion 

It is not surprising that the astronomical PAH model is able to reproduce the vdB 133 UIR 
spectrum. It is also not unexpected that the PAH mid-IR emission spectra in regions devoid 
of UV radiation resemble those in UV-rich regions. This can be readily seen from Figure 7 of 
Draine & Li (2001) which shows that a 0.5/im photon will excite a Ac = 100 (mean size for the 
Milky Way mixture) PAH molecule to a vibrational temperature of T ~ 390 K, resulting in strong 
emission in the 6.2-11.3/im UIR bands. For smaller PAHs, softer photons are able to fulfill this 
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excitation task. Therefore, the detection of UIR bands in UV-poor regions does not contradict 
PAHs as the UIR band carrier. The lack of a precipitous drop for cool stars (T c g < 10000 K) 
in [XI \ (12 (im)] iras / ^boi (FIR) simply indicates that interstellar PAHs have absorption extending 
into the visible, presumably as the combined result of somewhat larger sizes than those currently 
studied in laboratory, plus modification of the absorption properties by ionization. 

Papoular (2000) proposed that PAH excitation could result from energy released during PAH- 
catalyzed H2 formation. While this excitation mechanism may be present, it does not appear able 
to provide sufficient excitation to account for the observed UIR emission (Pagani et al. 1999). In 
any case, our model calculations show that excitation by starlight can account for the observed 
UIR intensities. 

We note that our approach to estimating the absorption properties of astronomical PAHs was 
simplified (Li & Draine 2001) in that we did not distinguish negatively-charged anions (PAH - ) 
and multiply-charged cations (PAH" + ) and anions (PAH™ - ) from singly-charged cations (PAH + ). 
For example, laboratory studies of the electronic absorption spectra of coronene (C24H12) in so- 
lution have shown that its di-anion (C^Hj^) has markedly enhanced visual/near-IR continuum 
absorption than its mono-anion (C24H 12 ) and neutral counterpart (Hoijtink 1959). However, the 
fact that our simplified PAH model can successfully reproduce the vdB 133 UIR spectrum and 
the T e ff dependence of [AiA(12/im)]iRAs/^boi(FIR) suggests that the model provides a fairly good 
approximation to the real absorption properties of astronomical PAH mixtures. 

6. Conclusion 

We have modeled the excitation of PAH molecules in UV-poor regions. It is shown that the as- 
tronomical PAH model provides a satisfactory fit to the UIR spectrum of vdB 133, a reflection nebu- 
lae with the lowest ratio of UV to total radiation among reflection nebulae with UIR bands detected. 
It is also shown that astronomical PAHs can be pumped by cool stars with even less UV radiation. 
It is further shown that the PAH model predicts a dependence of [A/A(12/um)]iRAs/^boi(FIR) for 
reflection nebulae which is consistent with observations for 3000 K < T c g < 30000 K. We conclude 
that PAHs appear able to account for the UIR band emission observed in reflection nebulae. 

We thank K.I. Uchida for providing us with the ISO spectrum of vdB 133, F. Allard, C. 
Dominik and P.H. Hauschildt for their help in obtaining stellar model spectra, and R.H. Lupton 
for the availability of the SM plotting package. We thank L.J. Allamandola, J.M. Greenberg, K. 
Sellgren and K.I. Uchida for helpful discussions. This research was supported in part by NASA 
grant NAG5-10811 and NSF grant AST-9988126. 
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